UMN-TH-1906/00 
TPI-MINN-00/26 
GRP/00/541 

in pre-big bang string cosmology hep-th/0006054 



Reheating and dangerous relics 



Alessandra Buonanno a , Martin Lemoine 6 and Keith A. 01ive c 
a Theoretical Astrophysics and Relativity Group, 
California Institute of Technology, Pasadena, CA 91125, USA 
b DARC, UMR-8629 CNRS, Observatoire de Paris-Meudon, F-92195 Meudon, France 
c Theoretical Physics Institute, School of Physics and Astronomy, 
University of Minnesota Minneapolis, MN 55455, USA 
We discuss the mechanism of reheating in pre-big bang string cosmology and we calculate 
the amount of moduli and gravitinos produced gravitationally and in scattering processes 
of the thermal bath. We find that this abundance always exceeds the limits imposed by 
big-bang nucleosynthesis, and significant entropy production is required. The exact amount 
of entropy needed depends on the details of the high curvature phase between the dilaton- 
driven inflationary era and the radiation era. We show that the domination and decay of 
the zero-mode of a modulus field, which could well be the dilaton, or of axions, suffices to 
dilute moduli and gravitinos. In this context, baryogenesis can be accomodated in a simple 
way via the Affleck-Dine mechanism and in some cases the Affleck-Dine condensate could 
provide both the source of entropy and the baryon asymmetry. 



I. INTRODUCTION 

In the pre-big bang scenario (PBB) Q, the standard Friedman- Robertson- Walker (FRW) post-big bang 
picture emerges as the late-time history of a Universe which, in a prehistoric era (the so called pre-big bang 
era), underwent an inflationary expansion driven by the growth of the universal coupling of the theory. This 
latter phase is also referred to as dilaton driven inflation (DDI) as its super-inflationary dynamics are driven 
by the kinetic energy of the dilaton field. A crucial difference between the pre-big bang model and standard 
inflationary theories, is that in the PBB, the universe starts its evolution in a classical state, the most general 
perturbative solution of the tree-level low-energy string effective action. The analysis of this initial state, 
and its naturalness, has led to a debate, as to whether the initial conditions needed to solve the horizon 
and flatness problems can be deemed natural Q . The problem of the graceful exit of the inflationary era is 
also an unsolved question. No-go theorems preventing the branch change from DDI to the dual solution of 
the FRW type have been demonstrated when either an axion field and a dilaton/axion potential or stringy 
fluid sources have been introduced in the tree- level effective action [[|. It is now recognized that if the 
branch change from DDI to FRW expansion is to occur, it should arise as a consequence of quantum loops 
effects and/or high curvature corrections, and encouraging progress has been made in this direction [Q, j^j. 
Nevertheless, the PBB model remains an attractive variant to standard inflationary cosmology, notably since 
the initial state of the Universe lies in the weakly coupled regime of string theory, and its dynamics are thus 
well controled by the tree-level low-energy effective action. This is in contrast with standard inflationary 
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theories which generically experience difficulties in extracting a well-suited Lagrangian for the inflaton field 
from a well-defined underlying fundamental theory, and in which the initial inflaton field values are usually 
of order the Planck scale. 

In recent years, significant effort has been spent on understanding the physics of the pre-big bang initial 
state and its high curvature phase, and on extracting observational predictions for this scenario. In this 
respect, one should note the prediction of a stochastic background of gravitational waves || , whose amplitude 
might be well above that predicted by models of standard inflation, as well as the amplification of quantum 
vacuum electromagnetic fluctuations due to the non-conformal coupling between the gravitational and 
the electromagnetic fields. More recently, it has been shown that the amplification of (universal) axion 
quantum fluctuations might provide adequate seeds for the formation of large-scale structures, and the 
resulting large and small angular scale anisotropies of the cosmic microwave background have been calculated 
. This provides a characteristic signal of non-Gaussian isocurvature perturbations^. So far little attention 
has been paid to the phenomenology of the post-big bang FRW era, and notably on the mechanism of 
reheating. It has been proposed that reheating could proceed via gravitational particle production |Tl| . 
However, as we argue here, this predicts the presence of too many dangerous relics, much like non-oscillatory 
inflationary models ]l2| ], and notably an abundance of gravitationally interacting scalars (e.g. moduli) well 
in excess of the limits imposed by big bang nucleosynthesis (BBN) on the abundance of late decaying massive 
particles. On top of that one naively expects that in some PBB scenarios the FRW era starts at a high 
Hubble scale, of order the string scale ~ 10 17 GeV, and therefore gravitinos and moduli should also be 
produced copiously in scattering processes in the thermal bath. 

These considerations warrant the present detailed study of the physics and phenomelogical problems of 
the PBB scenario in the post-pre-big bang era. We will find that indeed, for the variants of the PBB scenario 
hitherto proposed, there is inevitably need for significant entropy production. However, as we will argue, 
there exist various possible and natural sources of entropy production in pre-big bang models, notably the 
domination and decay of the zero-mode of a modulus, of the dilaton, or of axions, depending on the masses of 
these fields, which can dilute effectively moduli, gravitinos and monopoles. We will also show that this allows 
one to efficiently implement Affleck-Dine baryogenesis. In the following, we will thus focus on reheating, on 
the gravitino/moduli problem, and on the origin of the baryon asymmetry of the Universe. We will try to 
remain as general as possible, in particular with respect to the possible presence of an intermediate phase 
between DDI and FRW. 

The outline of this paper is as follows. In Section II we review the dynamics of the PBB model, and 
some variants (as far as the intermediate phase is concerned) proposed in the literature. In Section III, 
we provide a book-keeping of the particle content of the Universe at the beginning of the FRW era, and 
explicitly show the need for entropy production. In Section IV we review the various possibilities for entropy 
production in the context of the PBB model when the transition from DDI to FRW occurs suddenly, and 
discuss baryogenesis. We defer the study of intermediate phases to Sections V and VI, since the consequences 



1 One should note here that the recent high precision small angular scale data of the BOOMERANG and MAXIMA 
experiments Jl(J do not seem to confirm the predictions made in Refs. 



in that case are different but the logic of the argument is the same. We summarize our results in Section 
VII. 



II. THE PRE-BIG BANG ERA 



Let us first start by reviewing the different eras and dynamics envisaged in the PBB model. We shall 
restrict ourselves to the four-dimensional tree-level low-energy string effective action derived from heterotic 
string theory compactified on a six torus |l3|JI^ ] , whose bosonic sector is described by: ^ 

S oS = ^| J d 4 x^g~ [K + gT d^ip d u if - gT 8, t a a d v a a + £ ma ttcr] , (2.1) 



where A s = \V = \/8tt/M s is the string-length parameter, M s denotes the string mass and the following 
relations hold: 



Henceforth, we shall assume that at the present time go ~ 0.01—0.1. We restrict the internal compact space 
to a diagonal metric g a b = e 2<T ° 5 a i, with a, b = 4, . . . , 9, and we denote by ip the effective four-dimensional 
dilaton field ip = <f>io — ^2 a a a- The matter Lagrangian, £ ma tter, is composed of scalars, gauge fields and 
axions. We assume that the gauge and axion fields do not contribute to the cosmological background, i.e. 
we deal only with their quantum vacuum fluctuations. For the gauge fields, we consider the heterotic gauge 
field and the Kaluza-Klein gauge fields related to internal components of the metric and the three-form 
H^p, respectively V* and Wf 

/Igaugc fields = \ a' F M „ F» v - i a' e 2 °« V^ a - X -ol e~ 2a ° W^ a W% , (2.3) 

where 

F ia , = d li A u -d v A ll , Vpf = dfj,V u a — d v Vp , W^ a = d„W va - d u W^ a . (2.4) 

Finally, for the axion fields we have 

A™ - ~\ e 2v gT d^Ad v A - 1 gT e" 2 ^ e" 2 '" d^B bc d v B hc , (2.5) 

where B a j, is the pseudo-scalar field associated to the compactified components of the anti-symmetric field 
living in ten dimensions, while A is the axion related to the anti-symmetric tensor H^ vp in four dimensions by 
the usual relation H^ vp = e A " yp<T e v d„A (where t^ vpa is the covariant full antisymmetric Levi-Civita tensor). 

Let us consider first the simplest scenario in which the universe undergoes a super-inflationary evolution 
up to conformal time rji, at Hubble scale H± ~ M s and where the radiation dominated era is supposed to 



2 We use the conventions ( — ,+,+,+) and TV* vpa = T% a >p — ... , 1Z pv = 1Z P flpu . Units are % = k = c = 1, 
mpi = A/ P1 /(87r) 1/2 ~ 2.4 x 10 18 GeV is the reduced Planck mass 
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start, i.e. where the branch change from DDI to FRW occurs. The cosmological background during such 
DDI era is given, in conformal time, by 

a( V ) ~ (- V ) 5 /^ , <p ~ Vl + ^—j log(-r?) , (2.6) 

<r a ~ <r a i + r^l hg(-ri) , 5<0, a >O, (2.7) 

(note that with respect to the cosmological time a(t) ~ (— t) s ). A scalar field x with canonical kinetic term 
evolves classically during DDI as: 

X = -£= tog(-Tj) , (2.8) 

where the parameters S, f3 a , and j3 x satisfy the Kasner-type constraint: 

l = 3<5 2 + ^/32 + i^. (2.9) 

a " 

This Kasner constraint can be rewritten as a relation between S and an effective set of a which parametrizes 
the evolution of other scalars, including internal moduli. We will thus neglect (3 X in the following. We have 
also assumed that the four-dimensional non-compact space-time expands isotropically, while the contraction 
of the six internal dimensions can be anisotropic. After the branch change has occurred, the metric is that 
of a spatially flat FRW space-time; at that point the kinetic energy of the dilaton has become negligible, 
and the dynamics are thus driven by radiation, so that a(r) > 771 ) oc 77. The ulterior evolution of the dilaton 
is an unsolved question. We will assume that the dilaton is fixed in the radiation era p5| , but we will also 
indicate explicitly the dependence on the string coupling 171 (corresponding to the value of the coupling at 
the start of the radiation dominated phase) in our results. In particular, in the radiation era, the critical 
energy density p c as a function of the Hubble scale H is: p c — ( 3/ 871- ) i?., - 2 i? 2 M 2 . 

In some pre-big bang scenarios, the branch change from DDI to FRW is not instantaneous, and one con- 
siders an intermediate phase whose dynamics are obtained by taking into account higher order corrections 
to the low-energy effective action, such as finite size string effects and quantum string-loop effects. Un- 
fortunately, a thorough knowledge of the dynamics and duration of this intermediate phase is still lacking. 
Cosmological solutions, which partially describe the high curvature phase, have nevertheless been proposed 
in the literature, most notably: (i) the "string" intermediate era, obtained by solving the equations of motion 
with only the first order corrections in a' included Q] , and (ii) the "dual-dilaton" intermediate phase [[l4| (see 
also |p^| where this scenario was discussed in the more general framework of non- minimal models ), where 
one assumes that all a' corrections are sufficient to provide by themselves (without including string-loop 
effects) a sudden branch-change from the DDI to another duality-related vacuum phase of the FRW type. 
In the "string" intermediate phase, the Hubble parameter is constant hence the dynamics is inflationary in 
the string frame, while in the "dual-dilaton" era, the Hubble parameter makes a bounce around its maximal 
value at the string mass. To simplify the discussion, we shall often assume in both cases that the internal 
dimensions have been stabilized in some way before the Universe enters the intermediate era. We fix at rj s 
the time and H s the Hubble scale at which the Universe transits from the DDI era to the intermediate phase. 

For the "string" intermediate phase, one obtains j|: 
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<fl)~—J-, tp(v)~<p s -XIog(- V ), C>0, (2.10) 
H s rj 

where £ is an arbitrary parameter which governs the growing of the dilaton field, while with the "dual-dilaton" 
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a{v) ~ (-V) e/(1 - e) , <P~V. + Y=}l°g(-ri), (2.H) 



and 8 satisfies a Kasner constraint similar to Eq. (12 



III. THE POST-PRE-BIG BANG ERA 

In this section, we will consider the simplest version of the pre-big bang scenario with a sudden branch 
change from DDI to FRW, i.e. no intermediate phase of dynamics. 

A. Particle content due to gravitational production 

The particles present at the very beginning of the radiation era result from gravitational particle pro- 
duction, in contrast to standard inflationary models, in which the post-inflationary era is dominated by 
inflaton condensates, which later decay into radiation in the reheating process. In fact, in the PBB scenario 
it is the kinetic energy of the dilaton, which drives the DDI phase, that is converted into gravitationally 
created particles, whose energy density will drive the FRW era. One can provide a simple estimate of the 
energy density contained in fields subject to gravitational particle creation, when no intermediate phase is 
present. If a\ = a(r]i) is the scale factor at the branch change, then r/i — (aiHi)^ 1 (since H = a'/a 2 ), 
and fci = 1/r/i represents a comoving wavenumber corresponding to the horizon size at the branch change. 
We also define dpj/dhnk as the energy density spectrum in particle species j as a function of wavenumber 
k. Then one obtains dpj/dhik ~ for wavenumbers k > k%, since those modes have remained within the 
horizon at all times, and could not be excited by the gravitational field. For fluctuations that exited the 
horizon during DDI and re-entered during FRW, i.e. those modes with wavenumber k < k\, one generically 
obtains dpj/d\nk oc {k/kx) n o , and "-j 301 is the spectral index acquired by species j due to the dynamics 
of the DDI phase and transition into FRW. One imposes n® 1 * 1 > so as to avoid infrared divergences, i.e. 
large-scale inhomogeneities (see also below) and the energy density in species j is dominated by the energy 
density in the log interval around k\, so that pj ~ (dpj/dhik)\ k . Moreover, k\ corresponds to the maximal 
amplified wavenumber: this mode has exited and re-entered the horizon at the same time, and roughly one 
particle has been produced in that mode. Gravitational particle production thus respects a democracy rule 
PJ, namely all species share roughly the same energy density ~ (ai/a) 4 iJ 4 , corresponding to one particle 
produced with momentum Hi in phase space volume ~ Hf. Therefore, for all species j, pj ~ (ai / a) 4 H\ at 
times rj > rji, and consequently flj ~ g 2 (Hi/M s ) 2 , where £2j denotes the density parameter in species j and 
gi = e Vl / 2 is the value of the string coupling at the beginning of the radiation era [Mpi(?ji) = M s /gi\. 

When more accurate calculations are performed, one finds that the above democracy rule is satisfied to 
within less than an order of magnitude between different species, and one obtains (with 7jP di ^ 0) 
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and 



3?r \M s 



fJi 



1 



v > m , 



(3.1) 



(3.2) 



where Afj is the number of helicity states in species j. In Eq. (3.2), fcu is the wavenumber corresponding to 
the horizon size at time r\, i.e. fcn = 1 / 77. In effect, only modes whose wavelength is smaller than the horizon 
size can be thought as propagating as particles, and can be included in the energy density. For modes whose 
wavelength is larger than the horizon size, the definition of an energy density becomes gauge-dependent. 
Nevertheless, since we impose ?tP di > to avoid infra-red divergence problems, the contribution from the 



term in the absolute value in Eq. (3.2) is negligible for r\ rji, and the density parameter reduces to that 
deduced above by heuristic arguments, up to the fudge factor 4A/}/3tt. At this point, one should note that 
some fields j, and notably the axion A, can actually have n® Dl < f\ (see Tab. |). Again, to avoid infra-red 
problems, we shall impose n^ 01 > 0, that is — 1/3 < <5 < [l9[| . For the particular case of the PBB dynamics, 
it has been shown that the fields subject to particle production are those of spin 0, 1 and 2. One should 
mention that in Einstein gravity, abelian gauge fields are conformally invariant, and thus not gravitationally 
amplified; here, their conformal invariance is broken by the time evolution of the string coupling. Fermions 
(spin 1/2 and 3/2) are not produced |^(| (see also Refs. |2l| Q), at least when effects of compactification 
are neglected (see below). 



Particles 



n (string phase) 



n (dual-dilaton phase) 



moduli 

axion A 
axion B a t 
Heterotic photons 
KK photons 



-0.46 -j- 1 
-l-=-3 
1^-3 
1^-3 




-2.9^-1.4 





0.54 
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1.28 



—0.73 -5- 0.54 



TABLE I. Spectral slopes for the particles which have been amplihed gravitationally during the PBB phase 
assuming non-dynamical internal dimensions during the intermediate phases. The spectral indices in the three 
columns refer to fluctuations that exited the horizon during DDI, during the "string" phase, or re-entered during the 
"dual-dilaton" phase, respectively. 



3 The fact that axion fields can have negative spectral slopes is not a prerogative of the heterotic string model under 
study, in fact Copeland et al. [jl7| have shown that in the type IIB string model, with three axion fields, one of them 
at least must have nf 01 < 0, which can pose serious problems for the PBB model. On the other hand it has been 
shown recently |l^] that with a SL(A, 5R)-invariant effective action, there exists a region of parameter space where all 
the axions have nf m > 0. 
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In Tab. | we summarize the values of the spectral slopes for all of the particles present in the model 
assuming, for simplicity, non-dynamical internal dimensions during the intermediate phases. The range 
of values of n DDI have been obtained varying 8 and considering the possibility of having either one or six 
internal dynamical dimensions during DDI. Note that the spectral slope n DDI for the moduli fields have been 
obtained while considering them as part of the background ||l9|j , while in the determination of n with an 
intermediate phase, we neglect their presence in the background. Spectral slopes for other scalars depend in 
principle on their kinetic terms, and for simplicity we will assume that they have the same slopes as moduli 
fields, i.e. that they have canonical kinetic terms. The spectral slopes n in the second and third columns 
correspond to the slopes for fluctuations that exited during the "string" phase and re-entered during FRW, 
or exited during DDI and re-entered during the "dual-dilaton" era, respectively, and will be discussed in 
Sections VI and V. 



B. Thermalization and reheating 

In this section we analyse the thermalization and reheating process due to gravitational particle production. 
Let us start by considering the simple generic case with no intermediate phase. At Hubble scale Hi , all fields 
are produced with similar energy density ("democracy rule"). Let us denote by fif r the number of degrees 
of freedom in spin and spin 1 fields charged under the gauge groups of the observable sector. Similarly, 
if A/tot denotes the total number of degrees of freedom in spin 0, 1 and 2, i.e. that of the fields produced 
gravitationally, then the democracy rule implies that the fraction of energy density contained in radiation 
(in the observable sector) is f2 r ~ M r /A/tot- If the number of particles charged under gauge group is much 
larger than the number of gauge singlets, we get Cl r ~ 1. However, in some string models, the number of 
gauge singlets may actually exceed the number of charged states, and in this case, one would generically 
expect 0.01 < O r < 1. Henceforth, to keep the discussion generic we shall explicit the dependence on fi r . 

All fields carry typical energy E ~ Hi(a\/a), and the radiation number density n r ~ Q r p c /E. Gauge 
non-singlets interact with cross-section a ~ or /E and thus, thermalization occurs when the interaction rate 
T r = n r av > H (v ~ 1 is the relative velocity), i.e. at scale factor ath- 



Qth 

ax 



12 if Hi 



1 , m- l gta , , . 

r »1 



(3.3) 



For fi r ~ 1, <7i ~ 10 _1 , a = g\j '(4-71") ~ 10~ 3 and H\ ~ M s , thermalization occurs in ~ 9 e-foldings of the 
scale factor at H t h ~ 10~ 8 Hi. Let us observe that the value of a we have used differs from chgut ~ 1/20; 
this discrepancy, which is linked to the difference between the string scale M s and the GUT scale, is usually 
attributed to threshold effects. Note also that the various gauge fields A^, V? and we introduced in 



the action, Eq. (2.3), can in principle have different gauge couplings depending on the compactification. For 
simplicity we assume a single coupling constant, a = gf/An, which refers to the field. 

Even before thermalization is achieved, one can define an effective entropy density s — (4/3)p r /T r , where 
T r = (SO/vr 2 ) 1 / 4 ^* 1 ^ 4 pl^ 4 is an effective temperature, with g+ the number of degrees of freedom in the 
radiation after thermalization, i.e. including spin 1/2 fields that were not produced gravitationally but were 
re-created in scattering processes. This effective entropy will reduce to the standard entropy of the radiation 
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once thermalization has been achieved, and s ~ 0.2gl^ 4 flr^ 4 (HiM s /gi) 3 ^ 2 (ai/a) 3 . Then, if fl r ~ 1, reheating 
is complete once radiation has thermalized. If O r < 1, reheating would only be achieved once the fields that 
carry the remainder of the energy density have decayed to radiation. Such processes are constrained by big 
bang nucleosynthesis, which requires that at temperatures T r < 1 MeV, Q r ~ 1 to within a few percent. 
Nevertheless, as we will argue in the following subsections, it will be necessary to release a vast amount of 
entropy to dilute the dangerous relics produced. This entropy production may be viewed as a period of 
secondary reheating. 

C. Dangerous relics 

Using the above results, one can determine the number density of scalar fields with gravitational interac- 
tions present at the beginning of the radiation era and analyse their possible phenomenological consequences 
on BBN. In what follows, we will denote such scalar fields generically as moduli. Moduli are produced 
gravitationally as argued above, and one also expects them to be produced in scatterings of the thermal 
bath at time r\ > r\\. We will inspect each of these effects in turn, and discuss moduli and gravitinos. 

1. Moduli 

We adopt the generic notation Yj = rij/s for the number-density tij to entropy-density s ratio of species 
j; the entropy density in radiation can be written as before s ~ g X J 4 f?J 4 ■ Using Eq. ( |3.1[ ) for p r , one 
can rewrite s ~ (2ir 2 )~ 3 / i gl/ 4 Mr ^ 4 H 3 {ai/ a) 3 . Note the dependence on Mr^ 4 which counts the number of 
degrees of freedom produced gravitationally, namely those of spin and 1 . Due to supersymmetry, obviously 
M r ~ <?*/2, since accounts for these latter and their supersymmetric partners. 

Similarly, the number density of moduli n m ~ p m /E m , with typical energy E m ~ H\{a\/a), hence from 



Eq. Q n m ~ {l/2n 2 )Hf( ai /a) 3 , and: 

Y« ~ 0.3,?- 1 , (3.4) 



where the superscript g on Y m refers to gravitational production. For g± ~ 10 2 — 10 3 , Eq. (3.4) gives 
Y^ ~ 10 -4 — 10 -3 , which is well above the bounds imposed by BBN on the abundance of late decaying 
massive particles see e.g., ^5|. Indeed, let us contrast these estimates with the upper limits on Y m imposed by 
big bang nucleosynthesis (BBN) from photon injection. When applied to the case of moduli and gravitinos 
whose lifetime ~ 7~ 1 Mp 1 /m 3 , where m ~ 0(100 GeV) denotes the modulus/gravitino mass, and 7 is a 
fudge factor for the decay width (7 ~ 1/4 in the case of the gravitino) these constraints become |26|j27[| : 
Y m < 10- 15 for m ~ 100 GeV, Y m < 10" 14 for m ~ 300 GeV, and Y m < 5 x 10" 13 for m ~ 1 TeV 0. These 
bounds assume that the gravitino/modulus decays into photons with a branching ratio unity. Results weaker 



4 Note that Holtmann et al. define Y m with respect to the photon number density n 7 , not s, and today s ~ 7.0n 7 ; 
also, the constraints quoted assume 7 = 1; for 7 / 1, the mass estimates apply to instead of m. 
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by ~ 1 — 2 orders of magnitude would be obtained if the modulus/gravitino decays only into neutrinos, 
since high energy neutrinos produce an electromagnetic shower by interacting with the cosmic neutrino 
background |2Sf| . Moreover, stringent constraints in the high mass range m > 1 TeV would also be obtained 
if hadronic decay is allowed |||]. Thus a safe and generic limit is Y m < 10 -13 , which corresponds to the 
celebrated limit on the reheating temperature Trh ^ 10 9 GeV in standard inflationary scenarios. When 
considering these limits and the above results for the PBB scenario, one realizes that entropy production to 
the level of at least ~ 8 — 10 orders of magnitude is required. 

Moduli are also created in scattering processes of the thermal bath. The total amount of moduli present, at 
times rj 3> T]i , can be obtained by solving the Boltzmann equation with adequate production and destruction 
terms, with Y m — YJ^ as initial condition at rj = r\\. This equation, when written as a function of radiation 
temperature T reads: 



dY s 
dT 



-\-b — ►m+c 



*m Y c (7c-\-ni—*a-\-b 



a , b 



S 

HT 



(3.5) 



where the superscript s on Y m refers to moduli produced by scattering processes. In the above equation, 
Y a ^b, c denotes the number density to entropy density ratio of species a, b, c, and a a +b^m+c is the cross-section 
of the process a + b — > m + c. Generically, a and b are relativistic, in which case Y a ^ ~ 0.3/ [note that the 
first sum in Eq. (3.5) is over degrees of freedom of a and 6]. Since we found previously Y 9 n ~ 0.3/g* ~ Y a ^, c , 
which corresponds to equilibrium with radiation, the Boltzmann equation implies dY/^/dT ~ 0, i.e. the 
production/destruction of moduli in scattering processes is negligible as compared to Y^, and the final 
Y ~ Y 9 

1 rn. 1 <m - 



2. Gravitinos 

It has been argued recently (20| that gravitinos should not be produced gravitationally in the PBB scenario, 
if the gravitino is effectively massless, i.e. if the superpotential (W) ~ in the DDI and FRW eras. 
During the DDI phase, one indeed expects (W) = in a simple model. However, compactification of 
internal dimensions during DDI or non-perturbative effects to stabilize the dilaton in FRW should lead to the 
appearance of a superpotential, which would break the above condition, and result in gravitino production. 
Unfortunately, the magnitude of this mass term is very model-dependent and one cannot really determine the 
amount of gravitinos produced gravitationally. However, it should be noted that if one gravitino is produced 
per mode around the branch change frequency, corresponding to saturation of Fermi-Dirac statistics, one 
would find Y^ 2 ~ 0.3/ g* per helicity state as in the case of moduli. 

In any case, gravitinos are produced in scatterings of the thermal bath, in the same fashion as moduli and 



the Boltzmann equation (3.5) can be used substituting m — > 3/2 etc. If Yy 2 ~ 0.3/ g+ per helicity state, 
corresponding to equilibrium, then as before dY^ji/dT ~ 0. However if, as advocated in Ref. [l8| , gravitinos 
are not produced gravitationally in the PBB scenario, then Yy 2 = 0, and an estimate of Y3/2 is given by 
integrating the Boltzmann equation, neglecting annihilation and co-annihilation channels. This neglect is 
justified as long as the final value Y" 3 / 2 *C l/<?*, i.e. as long as equilibrium is not reached. Thus one obtains, 
using the Boltzmann equation for gravitinos, the simple result 30 - 32| 
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Y 3 s /2 ^^Y a (3.6) 

evaluated at the end of the PBB phase, where T is the gravitino production rate and Y a denotes the number 
density to entropy density ratio in species a. For the particle content of the MSSM, the gravitino total 



production cross-section is E to t ~ 250a/mp 1 |3^,Q (see also |34 for finite-temperature contribution to 
the gravitino production cross-section), and therefore, integration of the Boltzmann equation (assuming 
radiation domination a oc 77) gives: 

Yi.^lOOagf^gl'^ri^]' 2 . (3.7) 



For gi - 10 _1 , a = £?/(47r) ~ 10" 3 , n r ~ 1, ~ 10 2 - 10 3 and H 1 ~ M s , one thus finds Y 3 s /2 ~ 
1(T 7 - 10~ 5 ; this justifies our neglect of the annihilation channels in the Boltzmann equation. In any case, 
the destruction terms would ensure that Y 3 / 2 would never exceed its equilibrium value, so that Y3/2 ~ 
min Yy 2 + Y^ 2 , 1/g* is a good approximation to the final abundance of gravitinos, independently of the 
value of Y£, 2 . 



It should be noted that Eq. (3/7) evaluates the number of gravitinos produced before radiation has ther- 
malized [see Eq. (|3.3| )1, at the Hubble scale H±. Moreover, if fermions are not produced gravitationally, then 
the only charged non-singlets present at scale Hi are those of spin and 1, and the gravitino production 
cross-section should be smaller, since only channels a + b — > ^3/2 + c involving a, b of spin or 1 should con- 
tribute. However, we do not expect this uncertainty to exceed an order of magnitude [[30"|-|32|] . Furthermore, 
we used a ~ 10~ 3 as before, corresponding to gi ~ 0.1, and we neglected the running of a between M s and 
Mgut- However, it is easy to check that calculating Y3/2 with parameters corresponding to the GUT scale 
(Hubble scale i/cuT, coupling gi and ogut), one would obtain the same result as above. This is because 
the higher cross-section at the GUT scale S tot oc «gut, compensates for the smaller Hubble scale Hqvt- 
Overall, we estimate the uncertainty in the calculation of Y^ 2 to be < 1 order of magnitude, and the final 
Y3/2 exceeds by far the bounds imposed by BBN, similarly to moduli. 



3. Monopoles 

Finally, it is important to mention that the PBB scenario also suffers from the usual monopole problem 
due to GUT symmetry breaking (see also Ref. |p7|). Assuming that p monopoles form per horizon volume 
~ (4:Tt/3)Hq 3 at GUT symmetry breaking, one finds that the density parameter in monopoles today is 

where T c - 0.5ftr /4 # c 1/S 'mp( 2 is the critical temperature of the phase transition, mja is the monopole mass, 
and h denotes the Hubble constant today in units of lOOkm/s/Mpc. Naively, one expects p ~ 1/8 by 
counting the number of field orientations per horizon volume that would give rise to monopoles. However if 
the radius of nucleated bubbles at coalescence is much smaller than the horizon volume, one could actually 
obtain p>l§. 
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IV. ENTROPY PRODUCTION AND BARYOGENESIS 



The previous section indicated the need for a major source of entropy production in PBB models without 
an intermediate phase of dynamics. This is a stringent requirement, but, as we discuss below, sufficient 
entropy can be produced to solve the moduli/gravitino/monopole problems. Furthermore, as we argue in 
Section 4.B, this provides a natural framework for implementing baryogenesis in the PBB scenario. 



A. Sources of entropy and dilution of dangerous relics 

The late decay of non-relativistic matter is a simple way to generate entropy. Consider in addition to the 
radiation background the presence of matter with an equation of state pi = wpi and w < 1/3. Let us denote 
the value of the scale factor at the time the energy density pi is equal to the radiation density, po , by adom 
corresponding to a Hubble scale -ffdom- For a > adom, the Universe will be dominated by p\ until its decay 
at <zrh corresponding to a Hubble scale Hrji- To show the explicit dependence on the scale factor, let us 
write po = .g*dom/5o/a 4 where g*dom denotes the number of degrees of freedom at ad om and p\ = p\ / a 3 ( w+1 \ 
Then at adom we have, g*dom/5o = Pia dom ^"' • Assuming instantaneous decay, we can denote the energy 
density of radiation produced in the decay by, P2 = g*RnP2/a where 5*rh denotes the number of degrees 
of freedom at reheating and 5*rh/02 = Pi a tai < ' W+1 ' > ■ 

If we call so, the entropy density contained in po at a^n, then sq = |ff*domPo^ 4 /a RH . Similarly, 
the entropy in the radiation produced by the decay is si = |ff*RHP2^ 4 / a RH = fs*RHPi^ 4a RH w+1 ^ 4 — 
^9*im9%tmPo^ 4 a im W+1 ^ 4 / a ^dom^ W+1 ^^ ■ ^ we assume that the entropy release is large, we can write 

/ n \J/„ \ [3-9(u>+l)/4] 
A S =£i^^£i=(-£i5H ) (4.1) 
SO So V^dom/ Vadom/ 

We can also express the entropy change in terms of the Hubble parameter using (-ffRH/-ffdom) 2 = 
{a Aom /a K nf {w+l) so that 

A S =(^H) f^i) 0-3/4. (4.2) 

\<7*dom / V -«RH / 

Note that we included explicitly in Sdom a factor of fll^ 4 , which accounts for the fact that the FRW era may 
be driven by relativistic fields, but not by radiation (meaning gauge fields of the observable sector) . We also 
assume that the dilaton is fixed to its present value, at the latest by the time of domination. 

Depending on the equation of state, the exponent [4 — 3(1 + w)]/[2(l + w)] takes values from 1/2 for w = 
(non-relativistic matter) to — * +oo for w — > — 1 (cosmological constant), which is what effectively happens 
in standard inflation. Entropy can also be produced in first order phase transitions, albeit to a modest level, 



generally not more than ~ 1 order of magnitude |46 



In the following, we will be interested in the case of domination and decay of oscillations of a classical scalar 
field x in its potential V(x)- In the present scenario we will assume that initially the field x 1S displaced from 
its low-energy minimum by an amount xo ~ -^S ~ m Pi- This assumption is reasonable so long as the energy 
scales we are considering are much larger than the mass of the scalar field. If we stick with canonical kinetic 
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terms for the moduli during the PBB phase and we appeal to no-scale supergravity models to describe the 
particle content at the beginning of the radiation era, the flat directions corresponding to the moduli are 
still preserved, at least at tree-level |47|] . When supersymmetry breaking occurs the moduli will get a mass 
and we assume that the potential takes the simple form V(x) — rn x X 2 /2. 

The dynamics of a scalar field in its potential in the expanding Universe are well-known: the field is 
overdamped, and remains frozen to its initial value xo as long as H > m x . For H < m x , the field oscillates 
with an amplitude oc a~ 3 / 2 . Provided xo ?S m Ph the field comes to dominate the energy density after having 
started oscillating; if, as before, domination occurs at Hubble scale H^ om , then for H^ om < H < m x , the 
amplitude of since the Universe is still radiation dominated; for H < -ffdom, its amplitude oc H. 

If we denote by a x the value of the scale factor when oscillations begin, then the oscillations dominate at 
Odom = ( m P\/xo) 2 0'x- The field decays when H ~ T x , where T x is the decay width of x; as before, we write 
Hrh = r x . Assuming that x has gravitational interactions, T x ~ r y x m x /Mp- [ , where 7 X is a fudge factor, we 
find that x' s decay at orh = 7x 2 ^ {xn / m?\) 2 ^ (m x / M-pij^ 4 / 3 a x . Inserting these expressions for ad om and 
cirh into Eq. (4.1) with w — 0, we get 

where we have set g*dom = 200, and 3*rh = 10. 

In the above, we chose to select a gravitational decay timescale for the x field, as it represents the most 
efficient source of entropy, and x 1S therefore the coherent mode of a hidden sector scalar or modulus. 
In principle it is possible to obtain more entropy production if m x < 10 6 GeV. However the reheating 
temperature, given by 

-1/4 / m„ \3/2 



T xRH ~0.6GeV 



1/2 (9*BB\-^ f ™ X y" (44) 

lx \ 10 ) Vl0 6 GevJ ' 1 ' 



should not be lower than ~ 10 MeV for BBN to proceed unaffected, which requires m x > 6 x 
10 4 7 x 1 ^ 3 GeV pH]. Furthermore gravitinos are re-created in x decay to the level of: Y3/2 ~ 
10" 13 (m 3/2 /lTcV) 2 (m x /10 6 GcV)" 3 / 2 , so that one should impose m x > 10 5 - 10 6 GeV [||. Finally if 
R— parity holds, one needs to achieve T x rh <; 1 GeV for annihilations of LSPs to take place efficiently 
enough to reduce its abundance to cosmologically acceptable levels Q. Overall, it seems that As ~ 10 12 
represents, within an order of magnitude, the largest entropy production that is compatible with cosmological 
bounds for a displaced oscillating modulus. 

If at time 771, Y m ~ 0.33" 1 , the final abundance of moduli is given by: 



Yrr. 



where we assumed, for simplicity, g+ ~ 200 for the number of degrees of freedom in the radiation bath at 
time rji , and O r represents as before the fraction of energy density stored in particles charged under gauge 
groups. For higher g±, the numerical prefactor is further reduced as (<7^/200) _1 . We found in the previous 
section that the abundance of gravitinos produced gravitationally and in scattering processes does not exceed 
the abundance of moduli, and therefore the above estimate provides an upper limit to the final abundance 
of gravitinos. As for the monopole density parameter today, it is given by: 
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(A) My (to^v) gs)~>- <") 

This result was already obtained in Ref. |5?J , which studied the dilution of the monopole abundance through 
moduli oscillations and decay. 

Therefore, an initial value xo ~ mp\ appears sufficient to solve the moduli/gravitino problem of the PBB 
scenario, and marginally sufficient with regards to the monopole problem. If xo ~ 10 17 GeV, it is still 
possible to dilute the moduli and gravitinos down to acceptable levels Y m ~ 10~ 13 , but not monopoles. It 
should be pointed out, however, that the number p of monopoles produced per horizon volume is uncertain, 
and furthermore, that annihilations of monopoles with antimonopoles have been neglected in the above 
calculations. As a matter of fact, in various patterns of symmetry breaking, it appears that the monopoles 
are tied by cosmic strings, in which case annihilation of monopoles and antimonopoles would be highly 
efficient, leading to a scaling regime with one monopole per horizon volume at all times [BSpjJ. In this latter 
case, there would be no monopole problem at all. 

It should also be noted that we did not mention the possible moduli problem associated with the co- 
herent mode of those moduli whose mass < lOTeV, even though we considered entropy production due 
to one such coherent mode with mass ~ 10 6 GeV. As is well-known, an initial displacement of order 
mpi from the low-energy minimum of these moduli potentials would lead to a cosmological catastrophe: 
a reheating temperature < IkeV and an enormous post-BBN entropy production > 10 16 fjoj . Once 
the modulus starts oscillating, it behaves as a condensate of zero momentum particles with abundance 
n m /s ~ 10 7 (m/100 GeV) _1 / 2 (xo/ m Pi) 2 ^r^ 3 ^ 4 , where to, xo> an d s denote respectively the mass, ini- 
tial vev and radiation fraction of energy density when the modulus starts its oscillations, and the entropy 
contained in radiation. Hence, the above source of entropy cannot reduce sufficiently the abundance of 
these moduli, unless xo J$ 10 -4 TOpi. Unfortunately there is no well accepted reason why at high energy, 
i.e. after the branch change, these moduli should lie close to their low-energy minimum, and this problem 
affects all cosmological models, not only the PBB scenario. Nevertheless, if the string vacuum is a point 
of enhanced symmetry, one would indeed expect the moduli to lie close to their low-energy minima at high 
energy scales J4l| . In this respect, there is an interesting difference between the PBB scenario and the stan- 
dard inflationary models. In effect, in this latter class of models, even if the coherent mode is not displaced 
at the classical level, the generation of quantum fluctuations on large wavelengths will generate an effective 
zero-mode on the scale of the horizon at the end of inflation, displaced from its low-energy minimum: the 
so-called quantum version of the moduli problem |^,^|. In the PBB scenario, since the spectrum of the 
fluctuations of a scalar field is very steep at very large wavelengths (spectral slope = 3), its amplitude is 
small enough not to regenerate, at a quantum level, the zero-mode moduli problem. 

Enhanced symmetry on the ground state manifold would not apply to the dilaton fil|] , and in this respect 
one could wonder whether the dilaton could not play the role of the x field above, while other moduli would 
be fixed to their low-energy minimum for the above reason. Since the dilaton drives the DDI phase with 
its kinetic energy, and since it is not expected to lie exactly at the minimum of its low-energy potential at 
the end of DDI, this possibility seems rather natural in the framework of the PBB scenario. Furthermore, 
it should be noted that indeed, in some realizations of gaugino condensation, the dilaton acquires a mass as 
high as - 10 6 GeV @. 
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Finally, another solution to the cosmological moduli problem involves thermal inflation (J44J, or more 
generally a secondary short stage of inflation at a low scale. Indeed if this period of inflation takes place at 
a scale H <C TO3/2, the effective potential of the modulus during inflation will correspond to its low-energy 
potential (i.e in the vacuum of broken supersymmetry), and the modulus will be attracted exponentially fast 
to its minimum. 



B. Baryogenesis 

The above source of entropy comes with a bonus, namely baryogenesis can be implemented in a natural way 
via the Affleck-Dine (AD) mechanism jf5|. As already discussed above, the string model we are implementing 
has many flat directions. Generically, in these vacua, the scalar quarks and leptons have non-zero expectation 
values and can be associated with a baryon number and CP violating operators. Supersymmetry breaking 
lifts the flat directions providing a mass to the condensate made of squarks and sleptons, the so called AD 
condensate. When the expansion rate of the Universe is of the order of the mass of the AD field, this field 
starts to oscillate coherently along the flat directions carrying the baryon number. Finally, the subsequent 
decay of the AD condensate generates the baryon asymmetry. 

As it will be useful for the subsequent discussion, we will briefly outline how this mechanism works. Let us 
denote by $ the AD condensate, and m$ its mass. We assume that $ is initially displaced by an amount $0 
from its low-energy potential, for reasons similar to those previously discussed. The energy density stored in 
the condensate is simply p$ = m\<f> . $'s begin to oscillate when H ~ m$ at a = a$ and come to dominate 
the expansion when a = cid om = (mpi/$o) 2fl *- After oscillations begin the amplitude of the oscillations 
decreases as a -3 / 2 . The decay width of the condensate can be written as T$ ~ 7$m| ) / < i> 2 , where $ is the 
time-dependent amplitude of <&, and $ <C mj>\ at decay; 7$ ~ a 2 /An is a fudge factor. Thus the condensates 
decay when a — a<2* = ($Q^ 3 /mp{%|/ 9 m$' 9 )a$. This is true so long as the Universe is dominated by $ 
oscillations at the time of their decay, thus requiring that $0 > m^rrip® . 

The baryon number stored in the condensate oscillations is given by 

, 3 



n B ^e\^m<s, 1 $q $ 2 = eA$ m$ 1 $q (— ) , (4.7) 



where A$ is an effective B-violating quartic coupling. The entropy produced subsequent to the decay of the 
AD condensates is roughly s ~ <7+hwP$ > so that the produced baryon-to-entropy ratio is 

\ *5/2 / \ 3/4 , 2 / \ 1/6 

n B -1/4 eAj> $ / f a$_\ -r,- x l±.\ v ( 7$ TOpl 



<i> 

where e is a CP-violating phase. As entropy is produced in the AD condensate decay, the abundance of 



moduli and gravitinos produced at the end of the PBB era will be diluted. Using Eq. (4.1) above, one can 
easily show that As is given by 

A S =(^) 1/ %,- 1/6 ^/ 4 ^|V' ^ 
V y*dom / m $ m pl 

and is only O(10 6 ) for m$ ~ 100 GeV; the corresponding reheating temperature is 
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TmiMVg: ^r . (4.10, 



Thus we see that unfortunately, the entropy produced in the decay of the squark-slepton condensate cannot 
by itself produce the required source of entropy. 



As one can easily see from Eq. (4.8), the Affleck-Dine scenario of baryogenesis tends to produce too large 
a baryon asymmetry, with tib/s ~ O(l) if the AD condensate dominates the evolution when it decays. One 
possibility to reduce the baryon asymmetry that is of interest in the present context, is the late entropy 
production, as pointed out in |46|-[48|]. Moreover, it turns out that in the cosmological scenario envisaged 
here, the Affleck-Dine scenario seems to be the only model of baryogenesis capable of producing the required 
baryon asymmetry. Indeed, since As > 10 10 is required, and since BBN indicates a baryon asymmetry 
Ub/s ~ 4— 7x 10~ n , if baryogenesis takes place before entropy production, one needs to achieve ub/s ~ 0(1) 
initially, and only the Affleck-Dine mechanism seems capable of such a feat. 

Let us now consider the combined effect of an AD condensate and the late decay of a moduli field. Since 
777$ ~ 0(100 GeV), the x ne ld above will start oscillating before <f>. In the PBB scenario with no intermediate 
phase, the value of the Hubble parameter at the end of the PBB phase is Hi ~ M s at a = a\ and it is much 
larger than the value of H when x would start to oscillate, which happens at H ~ m x ~ 10 6 GeV. Hence, 
before x starts to oscillate the Universe is in a radiation dominated era. 

As before, one can determine the epoch of x domination (at scale factor a^om) by setting p r = p x , 
using (a x /ai) 2 = Hi/m x (a x is the value of the scale factor when x' s begin to oscillate), giving: 
ai/idom = m l/ 2 Xo/^-i^ 2 ' m p\- An analogous equation can be written for $ and we see that provided 
777$ < 10 6 GeV(77i x /10 6 GeV)(xo/^o) 4 , X wm dominate the energy density before $ (afom > °dom) anc ^ 
this condition is satisfied for most values of the parameters we are interested in. Moreover $ decays before 
X- Therefore, since x dominates the evolution before and since x decays after $, the total entropy pro- 



duced remains the same as in Eq. ( hS ) above and is sufficient for solving the moduli and gravitino problems 
of the PBB scenario. 

We will assume that the x ne ld dominates the energy density before <& begins to oscillate at a$; we 
will relax this assumption further below. This condition is true as long as m$ < to x (xo/wpi) 4 , which for 
777$ ~ 10 2 GeV and m x ~ 10 6 GeV, gives xo > 10 _1 mpi, which is quite reasonable in our context. Thus we 
can relate a x and a$ through a x /a$ = (m<f>/m x ) 2 / 3 (xQ/'mp\)~ 2 / 3 - In this case, we can rewrite the baryon 
number stored in the condensate oscillations in terms of a x as 

77 B ~ £ A$ 777$^ ^(-Y- ^ H 3 , (4.H) 

V a ) mjmpj V a ) 

As before, we will consider the gravitational decay of x to proceed with a rate T x — r y x m x /Mp V %'s decay 
at cipji when H = T x or when (a x / a,Rn) 3 = l x m x Xo 2m p\ /(^ 7r ) 2 ■ Subsequent to decay, the Universe reheats 
to (7*rhTp[ h = /9 x ((7rh) = ^ x Xo( a x/ a RH) 3 and the entropy density is just (4/3)p x /Tr,h- Thus the baryon 
to entropy ratio is easily determined to be |47]] 

^4 3/2 

nB ~m,\ -.1/2 .-1/4 ® mx r419 v 

* 777| 777 p j 

For an effective quartic coupling A$ ~ m% / ($0 + Mx ) > corresponding to superheavy gaugino exchange of 
mass Mx, one finds: 
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n B _ m -A 1/2 ( ™* \-l \ / m x ^3/2 $g 

T" 3 - Xl ° £7 - llOOG^vJ fcj lioTGevJ WTW (413) 



3/2 

+ J ".v 



and we assumed g*RH = 10 for simplicity. 

If the $ field starts to oscillate before the \ &e\d dominates, i.e. if xo < 10 _1 mpi, then the r.h.s. of 



Eq. (4.13) should be multiplied by (m^/m x ) 1 ^ 2 (xo/ m P\) 2 - The baryon to entropy ratio in this case is: 



_^ 2X 10 £7x / (___)(___) ^j^. (4-14) 



Both values of u-q/s obtained, Eqs. ( 4.13 ) and ( 4.14 ), are too large, but as discussed in Ref. |47|j , it is likely 
to be reduced through various mechanisms. For instance, if baryon number violation is Planck suppressed, 
then A$ ~ m^/mpj, which would reduce the above asymmetry if $o ^ mpi- Furthermore, under certain 
conditions, non-renormalizable interactions may reduce the initial vev of the condensate down to possibly 
Mx ~ 10 16 GeV, and the baryon asymmetry would be reduced by ~ 2 — 3 orders of magnitude [4^| . Finally 
sphaleron processing of the baryon asymmetry in the electroweak phase transition can also lead to reduction 
of Ub/s, by as much as ~ 6 orders of magnitude. Therefore, it seems reasonable to conclude that a baryon 
asymmetry of the right order of magnitude can be produced in the Affleck-Dine mechanism in this scenario. 



Let us now discuss the implications of the presence of an intermediate phase between DDI and FRW on 
the above conclusions. As we shall see, one of the main features is that the "democracy rule" of gravitational 
particle production no longer necessarily applies in the presence of an intermediate phase, and the distribution 
of energy density among the various components may be drastically altered. In some cases, this will imply 
that gravitinos and moduli can be more easily diluted in entropy production. 



V. DUAL-DILATON INTERMEDIATE PHASE 



The dual-dilaton era is characterized by two wavenumbers ki and k s that correspond to the horizon 
size at the branch change between the "dual-dilaton" intermediate phase and FRW, and between DDI and 
the "dual-dilaton" era, respectively. During the intermediate phase (IP), modes re-enter the horizon [since 
wavelengths oc a(i]) do not increase as fast as the horizon size H^ 1 = a 2 /a'], and therefore k\ < k s . As a 
consequence, one still expects to produce roughly one particle per mode at the highest wavenumber k s , since 
it exited and re-entered simultaneously at time r] s . We will also assume H s ~ M s , since H s correspond to 
the maximal value of the Hubble scale. At times r\ > rji, one finds an expression similar to Eq. ( |3.l|) for the 
energy density of species j: 




where ny is the spectral slope for fluctuations that exited the horizon during DDI and reentered during 
the "dual-dilaton" intermediate era. We again imposed a low wavenumber cut-off fcn corresponding to the 
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horizon size, but it will not play a role for rj > 771, since > as before for all fields. The integrated 

energy density can be written as: 



Pj (v) 



3 TT 

2^" 



V > Vi 



(5.3) 



A. Particle content and reheating 



From Eq. (5.3) we obtain that at times rj > r/i, the energy density is dominated by the field with the most 
negative spectral slope ro* p , and the democracy rule does not apply, unless all n* p > 0. The spectral slopes 
of the various fields depend in a non-trivial way on the dynamics of the internal dimensions during DDI 
and during the "dual-dilaton" intermediate era (see Tab. ||). Since details can be found in Ref. JL4|, we will 
simply restate the relevant results, but extend the analysis to the case with dynamical internal dimensions 
during the IP, i.e. 6 ^ I/a/3. 
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FIG. 1. The contours delimit the values of the smallest spectral slope mnx,(n* p ) as a function of the parameters 



S and 6, defined in Eqs. (2.8), (2.1C). The dashed line separates the regions in which either the model-independent 



axion A or the axion B a b dominates. 



In Fig. [l], we show a contour plot of the smallest spectral slope minj(nj p ) in the parameter space 6 — 9; 
we recall that 8 and 9 characterize the evolution of external dimensions during DDI and the IP respectively 



[see Eqs. (2.8), ( 2. 1C )] . The dashed line separates the areas which refers to the model-independent axion 
A, and to the axion Bat, associated to the internal components of B^ v . The spectral slope of these latter 
depend sensitively on the evolution of the 6 internal dimensions during DDI and IP, and in Fig. |l|, we chose 
to show the case in which 2 internal dimensions compactify isotropically. In the following, we will examine 
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the various cases in which a dimensions compactify isotropically, and the other 6 — a are stabilized during 
DDI and IP. We find that the following two possiblities arise: 

(1) if \5\ ~ 1/3 and 9 is close to its maximal value, which corresponds to saturation of the Kasner 
constraint Eq. ( |2.9|) , i.e. stabilized internal dimensions, then the smallest slope is negative, and is carried 
by A. In this case, the model-independent axion carries all the energy density at rj > rji. However, when 
one assumes that compactification is anisotropic, i.e. a dimensions compactify isotropically and 6 — a are 
stabilized during both DDI and IP one obtains a slightly different picture. In some cases, as the one shown 
in Fig. ^, the minimal slope is everywhere negative, and the bulk of the energy density is carried by either 
A or B a t axions, depending on the value of 5 and 9. 

(2) if the smallest spectral slope is positive, meaning rt* p > for all j, then all fields roughly share the 



same energy density at time 77 > rj s . This case is therefore very similar to that envisaged in Section [II, 
i.e. for a sudden branch change. Notably, one finds Y m ~ 0.3/ g± as previously, and entropy production 
as before may be considered to eliminate this problem. We will thus ignore this case in the following, and 
rather concentrate on case (1), assuming that one of the axions dominate the energy density at time rji. 

In the following, we denote generically as A the axion field that carries the energy density at times rj > r]\, 
and when necessary we will specify whether A is A or B a \y. The fraction of energy density contained in 
species j is: 



m 




1 + 




-nf" 


Vs 




Vs 





?7 > Vi 



(5.4) 



where we used k\/k s = r\ s jr\\. One can derive a relation between the duration of the "dual-dilaton" phase 
\r]i/rj s \ and the coupling constant g\ at time 771 from the criticality condition p_4 ~ p c , which gives: 

(5.5) 



where the exponent e_4 = (29 — 1)/(1 — 9) + nj^ /2 and is negative in the region of the parameter space where 



X < 0. To derive Eq. (p.5[) we used a(rj) oc 1 77 



9/1-6 



H oc 1 77 



-1/1-6 



and g cx \rj 



(38-l)/2(l-6 



Eq. (O), and 



p c = (3/8n)g^ 2 H^M^ atr] = r] 1 . If A is the model-independent axion A, n 1 ^ = 4/(1 - S) - 49/(1 - 9), while 
if A is a B a b axion, n l J = 4 — 4\/3 a / (1 — S) — £ / (1 — 6)\, where f3 a and £ Q parametrize the evolution of the 
internal scale factor during DDI and the "dual-dilaton" phase, respectively; (3 a is tied to 5 by the Kasner 



constraint Eq. (2.9) and similarly for £ a as a function of 9 (see Ref. pM). Therefore Eq. (5.4) can be re- 



written as: flj ~ J^fj9x A [1 + .9i 



It is difficult to give quantitative estimates for the total fraction of 



energy density carried by radiation, since the various components of radiation have different spectral slopes, 
and some of them can be negative (see Tab. |). In the following, we thus assume that radiation can be 
considered, on average, as one species with number of degrees of freedom J\f r ~ <7*/2 (as before), with a 
positive spectral slope, which implies 



i7 r ~ J\f r <?"' 



(5.6) 



Note that this number can actually be of order 1 if Af r is sufficiently large as compared to g'l A , however, 



this case would be similar to a sudden branch change, since g x 



|?7iA?s| would be of order unity, and 



the results of previous sections apply. We thus assume fl r 1 in what follows. To be definite, let us take 
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gi ~ 10 1 , 9 = l/v3 (corresponding to stabilized internal dimensions during IP) we find, varying S in the 
range (0,-1/3): 

^~10- 7 -10~ 3 , Q r — N r (10~ 4 - 10 -3 ) , (5.7) 



while posing g\ ~ 10 



we get: 



lO" 20 - io- 8 , 



n r ~ Af r (io 



IO" 9 ), 



(5.8) 



In contrast to the scenario with no intermediate phase, when a dual-dilaton intermediate era is present, 
reheating cannot be provided by gauge non-singlets, because at r\ > r\\ they generically carry a small amount 
of energy (fij -C 1), as discussed above. Let us then investigate the possibility of reheating via the axion 
fields present in our PBB model. Reheating may proceed if the axion can recreate radiation by scattering or 
conversions with photons, that is through the processes A + A — > 7 + 7 and ,4 + 7 — > .4 + 7. The interaction 
rate of the latter channel is strongly suppressed relative to the the rate of the former, since the radiation 
number density is small. The interaction term between A and the gauge fields is of the form [A/M')!FJ- ', 
The mass scale M' ~ Mpi for the model- independent axion A, but M' can be different for the B a b axions, 
as it then depends on the compactification J49[] . Indeed, as shown in the action Eq. ( |2.5| ), the coupling of 
B a i, to WW or VV, depends on the expectation values of the internal moduli. 

The cross-section for A — A scattering thus is of the form a ~ E^/M' 4 , where the typical axion energy 
E A ~ Hi(ai/a), with H\ possibly much smaller than M s for a dual-dilaton intermediate phase. Finally, the 
axion energy density is given by: pj, ~ (3/87r)<7 1 ' 2 iJ 2 M 2 (ai/a) 4 , so that the interaction rate for scattering 
gives: 



Therefore scattering by the model-independent axion cannot provide reheating, as one normally expects 
Hi < M s . However, if (M'/M F {) < g{ /2 (Hx/Ms) 1 / 2 , then B ab - B ab axion scattering will produce radiation 
and reheat the Universe. 

If axion scattering is ineffective, reheating may still occur through axion decay, provided the axion mass 
is large enough to avoid problems associated with too low a reheating temperature. A typical axion lifetime 
is ta ~ M' 2 /m\, where vtla ~ A 3 /AfpiM' is the axion mass (neglecting fudge factors) ]49[| , and A denotes 
the highest scale at which gauge interactions to which A couples become strong. In particular, for A ~ 
10 14 GeV, corresponding to phenomenologically favoured scales of gaugino condensation in a hidden sector, 
m A ~ 10 6 (Af'/0.01M P1 )~ 1 GeV and t a - 3 x 10 -8 s (Af/0.01Af P i) 5 , and the reheating temperature is 
T RH - 10GeVg^ 4 (M70.01M P1 )- 5 / 2 . In the following subsections, we will discuss the implications of this 
on the moduli and gravitino problems. 



The estimate of the abundance Y m of moduli produced gravitationally can be obtained using the same 
methods as for the no-intermediate phase case. One actually finds the same result Y m ~ 0.3/(7*. This is 




(5.9) 



B. Dangerous relics 
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due to the fact that the spectral slope of the energy distribution of moduli and radiation is positive for all k 
momenta (see Tab. ||), and therefore the number density of moduli n m ~ (27r)~ 3 / 2 H^ 2 (a/ a s )~ 3 , where we 
used Eq. ( |5.l|) , and n m ~ p m /E m , with E m ~ H s (a/a s )^ 1 the typical moduli energy. One can also express 
the entropy density s ~ g X J 4 p^J 4 in a similar way, and obtain the above result for Y m . 

Just as in the no-intermediate phase case, one does not expect gravitinos to be produced gravitationally 
in the dual-dilaton phase if they are effectively massless. However, even if radiation has not thermalizcd, 
they can be produced by scattering during the dual-dilaton phase. Since the string coupling (hence the 
production cross-section) evolves with time during dual-dilaton phase, it is more convenient to write the 
Boltzmann equation in terms of conformal time, which when disregarding annihilations channels gives: 

dY s 

-^ = S«say a 2 , Vs<V, (5.10) 
arj 

where we recall that the scale factor a oc \rj\ e ^ l ~ d } during the IP, s oc a~ 3 , and S to t ~ 500<? 4 /Mg, and as 
before Y a ~ 0.3/ g± denotes the ratio of number density to entropy density of radiation quanta per helicity 
state. Using g 2 oc |jj|( 3fl_1 )/( 1_e ) during the IP, one easily obtains as before Y£, 2 — Y a T/H, with T ~ T,t ot Y a s 
the gravitino production rate, and the right hand side should be evaluated at time rji if 9 > 1/3, and at time 
rj s if 6 < 1/3. 

Thus, using (H tot s/H) m ~ 7Qg\ /2 gl ,A fir /4 (#i/M,) 1/2 , and Eq. Q, one finally obtains: 



with: 



^3/2-6^^ (^j , (5.11) 



Let us first observe that, if 9 < 1/3 the coupling g decreases in time during the IP (see Eq. (2.9)). In the 
scenario under investigation, i.e. with a dual-dilaton intermediate phase, it is assumed that quantum string- 
loop effects are never operative {g < 1) and that because of high-curvature corrections the Hubblc-parameter 
makes a bounce around its maximal value at the string scale. If we impose that the dilaton field reaches the 
present value at the end of the PBB era, unless we assume either an extremely low decreasing of the coupling 
during the IP or a very short intermediate phase (which would give for Y^ 2 roughly the same value as in the 
scenario with no-intermediate phase), we are forced to limit to the region of parameter space where 9 > 1/3. 



Having restricted ourselves to the case described by Eq. (5.12), we find that the exponent k takes values 
between ~ 2 and ~ 9 depending on the evolution of internal dimensions, and therefore, y 3 / 2 ~ 10 -4 — > 10~ 10 
for gi ~ 0.1, 5* ~ 100. For anisotropic compactification, in some region of the parameter space, one finds 
even higher values of k, for which there would be no gravitino problem at all. One thus finds that in the 
presence of a dual-dilaton phase, the gravitino is generically less efficiently produced (possibly much less) 
than in the no-intermediate phase case. This can be understood in the following way. If 9 > 1/3, the string 
coupling grows during IP, so that the cross-section S to t oc g 4 is very small at the beginning of the dual-dilaton 
phase, and gravitino production takes place at Hubble scale H±. However, unlike the no-intermediate phase 
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case, here one generically has Hi <C H s ~ M s , and therefore gravitino production is inefficient. Finally, the 
small number density of radiation quanta (recall H,.<1) also hampers gravitino production in this case. 



C. Sources of entropy and baryogenesis 

Due to the overproduction of moduli (and possibly gravitinos) , entropy production is still necessary in the 
present scenario. Entropy production is also necessary if the axion cannot reheat the Universe, e.g. if its 
coupling to the gauge fields is too weak, and its mass too small. However, there are several possible sources 
of entropy production, as we now discuss. 

If the axion field dominates the energy density at the beginning of the radiation era and it reheats 
the Universe through scattering or decay, reheating is accompanied by entropy production to the level of 
As ~ (p^/p r ) 3 / 4 , where pa and p r denote the energy density contained in the axion and in radiation at 
reheating. If the axion is still relativistic at that time, i.e. if it reheats the Universe through scattering at 
high energy scale, then As ~ Q r 3 ^ 4 , where tt r is defined as before at time r\i. The moduli abundance would 
be reduced to Y m ~ O.Sg^ 1 Clf and to obtain Y m < 10~ 13 , one needs Cl r < 10~ 13 . One can check that 
such a low value of fl r cannot be obtained for realistic values with g\ > 10~ 3 in the parameter space defined 
by 5 and 9, using the results so far obtained for the axion B a b (recall that the axion A cannot reheat by 
scattering). For example, in the case of the axion B a b with only 2 internal dimensions compactified and the 
other 4 stabilized, the smallest value of f2 r for g\ > 10 -3 is actually f2 r ~ 1CP 6 , with Af r ~ 100, which would 
imply Y m ~ 10~ 7 at the end of reheating and entropy production is still necessary. Moreover, note that 
tt r ~ 10~ 6 also corresponds to Hi ~ 10~ 13 M S , hence M'/M P i ~ 10" 8 , which is a rather strong requirement 
on M'. Nevertheless, one should recall at this stage that the domination and decay of an Affleck-Dine flat 
direction would lead to further entropy production As' ~ 10 6 [see Eq. (|4.9|) ]. Thus one could actually dilute 
the moduli and gravitinos down to acceptable levels with an Affleck-Dine condensate only (no \ modulus). 
This would however produce too large a baryon asymmetry, n^/ s ~ 0(1), unless the CP violation parameter 
is very small and/or electroweak baryon number erasure is very efficient |47|| . 

Let us now consider the alternate case, in which the axion acquires a large mass, and later decays to 
radiation. Entropy production then results from the decay of the axion oscillations around the minimum 
of its potential, whose generic form @ is V(A) ~ {A 6 /M^) [1 - cos(A/F)], with F ~ M'/tt 2 the axion 
decay constant. The axion vev Aq is in fact generically displaced by ~ M'/tt from its true minimum, and its 



coherent oscillations and decay will produce entropy. Using the results of Section III , one can easily obtain 



As ~ io r r^v (_^y um~V 4 . (5.i4) 

\M'J \0.01MpJ V10 14 GeVy r v ; 

The amount of entropy produced is only marginally sufficient, but it depends in a sensitive way on A and 
il r . For smaller values of the gaugino condensation scale, say A ~ few x 10 13 GeV, and, say f2 r < 10 -2 , 
reasonable values of As > 10 10 may be achieved. Note that baryogenesis can be implemented in the very 
same way as in Section III in this context. 

One should recognize that the above estimates remain somehow qualitative, since they make particular 
assumptions on the axion couplings, and therefore on the compactification process, whereas we assumed a 
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simple toroidal compactification. Nevertheless, our aim here is to show that there exists various possibilities 
to generate entropy to the level required. 

Even if axions cannot reheat the Universe through scattering or decay, one may still consider the mechanism 
discussed in Section III, where a modulus (the dilaton?) of mass ~ 10 6 GeV reheats through the decay of 
its coherent oscillations. The discussion is similar to that of Section III, up to the fact that there may be no 
radiation dominated era preceeding the x dominated era, if ffi < -ffdom (-ffdom denotes as before the Hubble 
scale at which x comes to dominate the energy density), i.e. if the dual-dilaton era ends as x dominates. 
However, as we now argue, this does not happen, and one always has Hi 3> -ffdom- Following the discussion of 
Section III to calculate -ffdom, using a cx H~ e in the IP era, one finds: -ffdom ~ m x[(^ n /^)(9iXo/M' s )] 2 ^ 2 ^ 36 \ 
where gi is the value of the string coupling at Hubble scale -ffdom, which marks the end of IP. In the absence 
of the x modulus, the transition to the radiation dominated era would take place at Hi ~ Hgg^ 1 ^^ 1 e%> 
[see Eq. ( |5.5[ )]. We consider that in either case, the value of the string coupling at the end of the IP, i.e. at 
ffdom or at ffi, should be close to its present value. One then must determine whether or not ffdom ^ ffi, 
and it can be checked that for nearly all values of the 5 — parameters, we have indeed -ffdom <C -ffi- The 
dual-dilaton phase thus ends before and radiation domination occurs before x dominates. All estimates 
made in Section III can thus be directly applied to the present case. It is interesting to note that here, one 
generically has f2 r <C 1 (the "radiation dominated" era is driven by A), and therefore entropy production 
is more efficient by a factor il^ 3 ^ 4 ~ 10 — 10 6 , for J\f r ~ 100 and gi ~ 0.001 — 0.1. The amount of moduli 
present at x decay is thus further reduced by this factor. The monopole abundance today is also reduced 
by a factor £l~ 1/4 - 2 - 100 [see Eq. §U§}. Baryo genesis can be implemented as before via the Affleck-Dine 
mechanism, and the baryon asymmetry is given by Eq. (4.14). 

In fact, if il r < 10~ 6 , entropy production to the level of As > 10 6 would be sufficient to dilute the moduli 
and gravitinos to acceptable levels, and such entropy could be provided by the Affleck-Dine condensate, in 
the absence of any x modulus. It should be noted, however, that such low values of fl r only arise when 
9i < 10~ 2 . 

As a conclusion, when either A or B a b axions dominate the energy density at the end of a dual-dilaton 
intermediate phase, there are various natural sources of entropy: axion scattering, axion decay, or domination 
and the decay of a modulus or an Affleck-Dine condensate. Since one generically has fl r <C 1, entropy 
production is more efficient than in the absence of an intermediate phase, and both moduli and monopoles 
can be diluted down to low levels. In some cases, the Affleck-Dine condensate provides enough entropy 
by itself to solve the moduli/gravitino problem, although one then has to cope with a very large baryon 
asymmetry from the decay of the condensate. 



VI. "STRING" INTERMEDIATE PHASE AND BLACK HOLES 

During the string era, modes exit and do not re-enter the horizon. Therefore k s < k±, and one expects to 
produce one particle per mode at ki since it exits the horizon and re-enters at the same absolute value of 
conformal time. 

However, the situation here is more delicate than in the case of the dual-dilaton phase. Indeed, a mode 
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that exits the horizon at conformal time r] ex , with rj cx < in the present scenario, will re-enter at time 
?7 rc ~ |77 ex | (if 77 ex <C 771 , i.e. if the mode exits well before the end of the "string" phase). This means that 
at time 771, which is supposed to mark the start of the FRW regime, only those modes with wavenumber 
~ fci, i.e. the highest frequencies, have re-entered. The modes that exited the horizon at the beginning of 
the "string" phase (wavenumber ~ k s ) will re-enter later, possibly much later, at conformal time |?7 S | 3> |?7i|. 
One can rewrite the energy density Eq. (3T) in this scenario, at time 77 with 771 < 77 < |t7 s |, which gives: 

Mi 



1 



771 < 77 < \n a \ , (6.1) 



and as before, we imposed a low wavenumber cut-off at the horizon size ~ I/77. In the case of the "string" 
phase, inspection of Table | reveals that the model-independent axion has the most negative slope = — 2( 
(C > 0), and its energy density pa ~ (AfA/2ir 2 )Hfaf\r]/T]i\ 2 ^ /a 4 . At time 771, only modes with wavenumber 
k\ have re-entered, so within the horizon, all fields share roughly the same energy density. However, as time 
goes beyond 771, the axion will quickly come to dominate the energy density. Assuming the dilaton field is 
fixed for 77 > 771, since H 2 oc pa, and H = a' /a 2 , it is straightforward to derive that a oc ?7 1+ ^. In this case, 
with £ > 0, the dynamics is driven by the axion fluctuations. This is inconsistent since the gravitational 
amplification of these fluctuations, in particular the spectral slope of the axion, were calculated assuming 
that the fluctuations would re-enter during a radiation dominated era, and with £ > 0, the expression for 
the scale factor shows that this is not the case. 

This inconsistency reflects the breakdown of the perturbative approach used to calculate the amplification 
of axion fluctuations. In effect, the calculation assumes that the quantum fluctuations can be treated as a 
perturbation on a fixed classical background, whereas in the present case, one should consider their back- 
reaction effect on the background spacetime. Moreover, since the axion field is assumed not to participate 
to the dynamics, its classical vev is zero, and the energy density stored in axion is pa ~ {d5A) 2 , where 5 A 
represents the fluctuation in the axion field. One should therefore include back-reaction up to second order 
in the fluctuations, in order to derive the dynamics of the era between times 771 and \ fj s \, and such an intricate 
calculation is well beyond the scope of the present paper. Note, however, that the above inconsistency does 
not arise for conformal times 77 > \rj s \, since nf Dl > for all j, and radiation domination should be a valid 
approximation. 

We thus consider, as an alternative, that black holes form on all scales comprised between k\ and k s . 
This is a possible outcome of the above dynamics, as black holes generically form copiously when relative 
overdensities of order unity re-enter the horizon. At late times 77 > \r) s \, the Universe will be dominated 
by those black holes that have not evaporated yet. The lifetime of a black hole Tbh oc Af 3 h oc H~ 3 : where 
Mbh denotes the black hole mass, and corresponds to the mass within the horizon at the time of formation 
at Hubble scale Hi. Therefore the Universe will be dominated by those black holes that formed last, 
i.e. on scale k s , and will reheat with the evaporation of those black holes. The Hubble scale H(\r) s \) ~ 
Hi |?7 s /?7i |~ 3 / 2 if the era between times 771 and \r] s \ is matter dominated (black hole domination). Then 
-^bh ~ 4irmpi[H (|77 s |)/?7ipi] (which corresponds to the mass within the horizon at that time) and the 
evaporation timescale of those black holes reads Tbh ~ nip^M^/mpi) 3 ~ (47r) 3 TOp 1 1 [i?(|77 s |)/77ipi] _3 , so that 
the black holes evaporate and reheat the Universe at a Hubble scale H cv ~ (47r)~ 3 i^ (|t7 s |) [i^ (|?7 s | ) / 7npi] 2 . 
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Such reheating by black hole evaporation was envisaged in Ref. [M in the context of the PBB scenario. 
Note that it is accompanied by entropy production, to the level of As ~ (pbh/ Pr-) 3 ^ 4 , where pbh and p r 
denote the energy densities contained in black holes and in radiation at the time of evaporation. This entropy 
production can be rewritten as As ~ 10 19 g\ /2 Q,7 3/4 [£T(|?7 s |)/10 7 GeV]~ 3 / 2 (iii/ M s) 1/2 , and the reheating 
temperature Trh ~ 0.6 GeV [H (\r] s \) / 10 7 GeV] 3 / 2 . Even though this entropy production may be sufficiently 
large to dilute the moduli created gravitationally during the DDI and "string" phases, the same moduli are 
also part of the Hawking radiation of the evaporating black holes. The number density to entropy density 
ratio of moduli and gravitinos present after evaporation in fact reads ]5l| : Y m ~ 10~ 9 [iJ(|77 s |)/10 7 GeV] 1 / 2 , 
assuming g+ ~ 200 at evaporation. Clearly, black hole evaporation eliminates one moduli problem, to 
reintroduce another, and further entropy production is necessary. 

Consider then a modulus x as introduced in Section IV. A, with mass m x . The amount of entropy 
produced by x 1S given in Eq. (|4.3|). However, if black holes have not evaporated by the time x would 



dominate (if black holes were absent), the r.h.s of Eq. (4.3) should be multiplied by {H cv / H^^) 1 / 2 , where 



-Hdom ~ w x (xo/ to pi) 4 corresponds to the Hubble scale at which % would dominate in the absence of black 
holes. Using H cv ~ 10 6 GeV[ff(|7/ a |)/2 • 10 15 GeV)] 3 , and assuming H cv < i?dom, the final moduli abundance 
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Y -lO- 17 ^ 2 ( m * V | H[lVsl) ) fi 3 / 4 (6 2) 

y m W 7 X ^ 10 6 GeV ; l v 2-10 1 5GeVy' T ' 1 ' 

hence, for Hi ~ 10 17 GeV, provided the phase in which black holes dominate does not last too long, i.e. 
\Vs/vi\ <4x 10 3 , then ff(|»7 a |) > 4 x 10 11 GeV and Y m < I0~ 13 . 

To summarize, black hole reheating of the PBB scenario, as envisaged in Ref. [po[ , does not solve the 
moduli problem; it also requires another source of entropy production, and the oscillations and decay of the 
X modulus would be sufficient, provided the "string" phase does not last more than ~ 8 e-folds of the scale 
factor. In this case, baryogenesis could also be implemented as in Section IV. B. 



VII. DISCUSSIONS AND CONCLUSIONS 

We find that pre-big bang cosmological models inevitably face a severe gravitino/moduli problem, as they 
predict a number density to entropy density ratio of gravitationally produced moduli at the beginning of 
the radiation era of the order of Y m ~ 0.3/g*, where g+ counts the number of degrees of freedom in the 
radiation bath at that time. These models also predict a similar amount of gravitinos, albeit somewhat 
smaller if gravitinos are not produced gravitationally during dilaton driven inflation, yet far in excess of the 
BBN bounds V m , 3 / 2 < 10~ 13 . 

Late entropy production, to the level of As > I0 5 — 10 10 depending on the details of the transition 
between the pre-big bang inflationary era and the radiation phase, is thus mandatory in the scenarios we 
have investigated. For the simplest pre-big bang model in which the transition is sudden, the amount of 
entropy needed is As > I0 10 and this is a strong requirement. However, sufficient entropy can be produced 
by the domination and decay of the zero- mode of a modulus field with mass ~ 10 6 GeV, initially displaced 
from the minimum of its potential by an amount M s . The Universe will start the FRW radiation era with a 
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temperature Trh ~ 1 GeV (see Eq. (4.4)). Moreover, the dilaton, which drives the pre-big bang dynamics, 
could also play the role of this modulus, as several scenarios of gaugino condensation predict a dilaton 
mass ~ 10 6 GeV, and since it can be generically displaced from its present value at the end of the pre-big 
bang inflation. Furthermore, this vast amount of entropy produced helps set the Affleck-Dine mechanism of 
baryogenesis in a natural framework, as it reduces efficiently the baryon asymmetry created in the decay of 
the baryon number carrying flat direction. Finally, it may also solve the usual monopole problem associated 
with GUT symmetry breaking, although this depends sensitively on the details of monopole formation at 
the GUT phase transition. 

We also examined variants of the pre-big bang model in which an intermediate phase of dynamics motivated 
by physics at high curvature takes place between the pre-big bang inflationary phase and the radiation era. 
In the case of the so-called dual-dilaton intermediate phase, one finds that the moduli/gravitino problem 
is still present, and entropy production is still necessary. However the problem of entropy production is 
relieved by the small fraction of energy density contained in radiation at the beginning of the radiation era. 
In effect, the energy density is generically contained in an axion field, either the model-independent axion 
A or internal axions B a b associated with the compactified components of B^. One can show that several 
natural sources of entropy may alleviate or solve the moduli/gravitino problem, notably the entropy produced 
in axion reheating via scattering (provided the axion decay constant F < 10 Afpi), or that produced in 
oscillations and decay of the zero-mode of the dominating axion, if the axion decay constant ~ 10~ 3 Mpi 
and its potential is generated by gaugino condensation in a hidden sector at scale A ~ few x 10 13 GeV. In 
some regions of parameter space (which parametrizes the evolution of internal and external dimensions), 
the amount of entropy needed to reduce the moduli/gravitino problem is sufficiently small (As ~ 10 6 ) that 
the domination and decay of an Affleck-Dine condensate can produce both the entropy and the baryon 
asymmetry of the Universe. In this case the hot big bang, which marks the beginning of the FRW radiation 
era, takes place at a temperature Trh ~ 10 5 GeV (see Eq. ( 4.10 )). 

In the case of the so-called "string" intermediate phase, one is at present unable to specify the dynamics 
of the era that follows the string phase (see Section VI). However, as we have argued, it is likely that 
microscopic black holes would form copiously. Black hole domination and decay produces entropy, which 
would dilute the moduli/gravitinos produced during the dilaton-driven and string phases, but moduli and 
gravitinos are also re-created in the Hawking radiation of the evaporating black holes. Here again, therefore, 
further entropy production is necessary, and the decay of a heavy modulus can produce sufficient entropy. 

At this stage we would like to comment on the implication of entropy production on the various predictions 
of the pre-big bang models. First of all the entropy production will not affect in any way the axion seeds 
of large scale structure considered by Durrer et al. In effect, as long as the axion perturbations lie 
outside the horizon, they are frozen, and do not suffer from the microphysical processes inside the horizon, 
i.e. they are not diluted by entropy production. More quantitatively, the density perturbation Spa/Pc in 
the axion field relative to the total energy density, can be written as a function of comoving wavenumber 
k and conformal time 77, as ||: &pa/Pc ~ (kri) 2 (Hi/Mpi) 2 (k/ki) n A D1 , for modes outside of the horizon, 
i.e. k < I/77, and where n^ 01 denotes the axion spectral index of density fluctuations. Clearly, when the 
mode re-enters the horizon, i.e. k — 1/r], for n 1 ^ 1 w (scale invariant spectrum), the density perturbation 
is independent of any entropy production. Note however that a nearly flat spectrum of axion fluctuations 
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corresponds to a region of parameter space in which the necessary amount of entropy release is large, of 
order As ~ 10 10 . 

Therefore, if perturbations of the \ held, which is the modulus responsible for entropy production, do not 
carry power on large scales, i.e. n^° l > 0, as would be the case if x were the dilaton for instance, the scenario 
envisaged by Durrer et al. Q for the axion seeds remains unaffected. In this framework, the pre-big bang 
predicts non-Gaussian isocurvature perturbations with a well defined signature in the small angular scale 
cosmic microwave background anisotropies. However, if the perturbations 5\ carry power on large scales, 
adiabatic perturbations would be produced on these scales at \ decay, since it dominates the evolution at 
that time. The study conducted in Ref. jl4j seems to indicate that the only fields in pre-big bang models 
that are liable to carry power on large scales are the axions A or B a b- In this context, the domination and 
decay of an axion, as considered in Section V.A., could lead to a novel scenario of generation of density 
perturbations in pre-big bang models. One should calculate carefully and examine the exact shape of the 
spectrum of metric fluctuations and their statistics, as it is known that axionic fluctuations are generally 
damped on large scales due to the periodic nature of the potential |5^| , 

Scalars generically carry steep blue fluctuations spectra in pre-big bang models, hence neither the mod- 
ulus x nor the Affleck-Dine condensate $ envisaged in Section IV are liable to produce long wavelength 
fluctuations at their decay. This is in some contrast to standard inflationary models, in which the decay of 
the Affleck-Dine field produces isocurvature long wavelengths fluctuations. 

On similar grounds, one does not expect that the spectrum of electromagnetic fields on the scale of the 
Galaxy should be diluted by entropy production, as these fluctuations were outside of the horizon at the time 
at which entropy was released. However, one expects that the relic gravitational wave background will be at 
least partly affected by the entropy release j53|,|o| . Whether this dilution affects the stochastic gravitational 
background in the range of frequencies which the upcoming experiment, LIGO, Virgo and LISA, are sensible 
deserves further investigation. 
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